MNRAS: in press 



Precessing collimated outflows in the 
planetary nebula IC 4846 

Luis F. Miranda 

Instituto de Astrofisica de Andalucia, CSIC, Apdo. Correos 3004, E-18080 Granada, Spain 

lfm@iaa.es 
Martin A. Guerrero 

Department of Astronomy, University of Illinois at Urbana- Champaign 1002 West Green Street, 

Urbana, IL 61801 USA 

marOastro . uiuc . edu 

Jose M. Torrelles 

Institut d'Estudis Espacials de Catalunya (IEEC/CSIC) and Instituto de Ciencias del Espacio 
(CSIC), Edifici Nexus, C/ Gran Capita 2-4, E-08034 Barcelona, Spain 

torrelles@ieec .fcr.es 
ABSTRACT 

We present [N n] and Ha images and high resolution long-slit spectra of the plan- 
etary nebula IC4846, which reveal, for the first time, its complex structure and the 
existence of collimated outflows. The object consists of a moderately elongated shell, 
two (and probably three) pairs of collimated bipolar outflows at different orientations, 
and an attached circular shell. One of the collimated pairs is constituted by two curved, 
extended filaments whose properties indicate a high velocity, bipolar precessing jet. A 
difference of ~ 10 kms -1 is found between the systemic velocity of the precessing jets 
and the centroid velocity of the nebula, as recently report for Hu2-1. We propose that 
this difference is due to orbital motion of the ejection source in a binary central star. 
The orbital separation of < 30 AU and period < 100 yr estimated for the binary are 
similar to those in Hu 2-1, linking the central stars of both planetary nebulae to interact- 
ing binaries. Extraordinary similarities also exist between IC4846 and the bewildering 
planetary nebula NGC 6543, suggesting a similar formation history for both objects. 

Subject headings: planetary nebulae: individual: (IC 4846) - ISM: kinematics and dy- 
namics - ISM: jets and outflows 
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1. Introduction 

IC4846 (PN G 027.6— 09.6) is a compact planetary nebula (PN) whose morphology has not 
been studied in detail yet. The only available information on its structure is provided by the VLA 
6cm continuum observations by Kwok (1985, see also Aaquist & Kwok 1990), showing several 
knots embedded in a faint elongated structure of ~ 3x2 arcsec 2 in size. The H/3 surface brightness 
(logF(H/3) = —11.34, Acker et al. 1992) suggests that IC4846 has a high electron density. This is 
corroborated by the small [S n]A6717,A6731 doublet ratio (Barker 1978; Acker et al. 1992) which 
reaches the limiting ratio for high electron density (iV c J;10 4 cm -3 ). Both its small angular size 
and high density indicate that IC 4846 is a young PN. The line intensity ratios (Acker et al. 1992) 
point to a moderate excitation ([N n]A6583/Ha ~ 0.07, [O ffl]A4959/H/3 ~ 3.7, and very faint He II 
A4686). IUE low-resolution UV spectra of the central star show a strong P-Cygni profile in the 
N v line 1 with a terminal velocity of ~ 1300 km s -1 . The occurrence of a strong stellar wind in the 
central stars of PNe is indicative of their early evolutionary stage (Cerruti-Sola & Perinotto 1985). 

The shaping processes of PNe occur early in their formation. Therefore, the signatures of 
these processes in the nebula are more recent and more likely to have withstood the fast nebular 
evolution in young PNe. This turns the observations of young PNe of chief importance in the un- 
derstanding of their formation mechanisms. Accordingly, IC4846 was included in our observational 
program of compact young PNe (see Miranda 1999 and references therein), and narrow-band Ha 
and [N n] images and long-slit high-resolution spectra were obtained to study its structure and 
internal kinematics. The results reported in this paper evidence the existence of bipolar, precessing 
collimated outflows in IC4846. 



2. Observations 

Direct images of IC4846 were obtained in 1997 July 24 with the Nordic Optical Telescope 
(NOT) at Roque de los Muchachos Observatory on La Palma 2 . The detector was a LORAL 15 
/im CCD at the HIRAC camera providing a plate scale of 0.11 arcsec pixel -1 . The filters were Ha 
(A ~ 6563 A, FWHM ~ 10 A), and [N II] (A ~ 6584 A, FWHM ~ 10 A). Exposure times 
were 900 s in each case. The spatial resolution is ~ 1 arcsec. The images were bias subtracted and 
flat-fielded following standard procedures within MIDAS package. 

Long-slit echelle spectra of IC 4846 were taken with the IACUB 3 spectrograph at the NOT in 
1 Based on INES data from the IUE satellite. 

2 The Nordic Optical Telescope is operated on the Island of La Palma by NOTSA in the Spanish Observatorio del 
Roque de los Muchachos of the Instituto de Astrofisica de Canarias 

3 The IACUB uncrossed echelle spectrograph was built in a collaboration between the IAC and the Queen's 
University of Belfast 
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1998 August 2. The detector was a Thompson 19 fim CCD with 1024x1024 pixels. A filter was 
used to isolate the Ha and [N n]A6583 emission lines in the 9 th echelle order. A Thorium- Argon 
spectrum was taken after each object spectrum for wavelength calibration purposes. Exposure time 
was 1800 s for each spectrum. The observed position angles (PAs) were 11°, 30°, 54°, and 140° 
with the slit centered on the object. The dispersion was 2.5 A mm -1 , and the achieved spectral 
resolution (FWHM) ~ 7.5 kms" 1 . The spatial resolution is ~ 1.5 arcsec. Absolute errors in the 
radial velocity are ± 1.5 kms" 1 , whereas relative errors (within each long-slit spectrum) are ± 0.5 
kms -1 . The data were reduced following standard procedures for long-slit spectroscopy within the 
IRAF and MIDAS packages. Thereafter the [N n]A6583 emission line will be referred to as the 
[N n] line. 



3. Results 

3.1. General structures 

Figures 1 and 2 present grey-scale and contour maps, respectively, of IC 4846 in the light of 
[N n] and Ha. The morphology of IC 4846, resolved in these images for the first time, displays a great 
wealth of structures. At the lowest intensity levels, the appearance of the nebula is elliptical with 
two bipolar features protruding at PA ~ 54°. These features, denoted A1-A2 (Fig. 2), are curved 
filaments enhanced in [N n]. A second, inner ellipse, observed in Ha and [N n] at intermediate 
intensity levels, can be recognized with its major axis at PA ~ 11°. Two compact knots, B1-B2, are 
observed at the tips of this inner ellipse. The innermost regions of IC 4846 are dominated in [N n] 
at the highest intensity levels by two bright structures, C1-C2, which are not clearly distinguished 
in the Ha image. An additional attached circular shell with a size of ~ 10" can also be identified 
in the Ha image. 

In order to emphasize the different structures of IC 4846 described before, we have obtained 
deconvolved images using the Lucy-Richardson algorithm as implemented in the MIDAS package. 
Deconvolved Ha and [N n] images to a resolution of ~ 0.5 arcsec (obtained with 30 iterations) 
are also presented in Fig. 2. These images clearly show the structural richness of the nebula. In 
particular, the curvature and filamentary appearance of A1-A2 can easily be recognized as well as 
the compactness of B1-B2. C1-C2 appear as point-symmetric structures elongated perpendicular 
to the radial direction from the center. In Ha the inner regions appear as an arc with an asymmetric 
intensity distribution. We note that the inner regions of the nebula resemble the morphology 
observed at 6 cm (Kwok 1985). In addition, the inner shell presents a very low ellipticity. 

Figure 3 shows position-velocity (PV) contour maps of the high resolution long-slit spectra 
of the [N n] emission line at PAs 11°, 54° and 140°. The spectrum at PA 30° is similar to that at 
PA 11° and is not shown here. The bipolar features A1-A2, B1-B2 and C1-C2 can be identified 
in the PV maps with point-symmetric pairs of components. We also note that no closed velocity 
ellipse can be identified in the PV maps of the [N n] emission line, which could be expected from 
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the elliptical structures observed in the images. The Ha spectra (not shown here) are dominated 
by emission from the inner regions. These regions show a tilt in the PV maps similar to the tilt of 
C1-C2 (Fig. 3). A1-A2 are clearly distinguished because of their high radial velocity. B1-B2 seem 
to be present, too, although the Ha line is not useful for a detailed kinematic study owing its large 
thermal width. 

A [N n]/Ha ratio map is also presented in Fig 1. A1-A2 and B1-B2 are highlighted due to 
their relatively strong [N n] emission. C1-C2 can also be recognized, although the line ratio seems 
to have lower values. A quantitative estimate of [N n]/Ha should be worked out using the long-slit 
spectra in order to single out both spatially and kinematically the different components. PV maps 
of the [N n]/Ha ratio for PAs 11° and 54° are shown in Fig. 3. A mean value of [N n]/Ha in the 
nebula of ~ 0.06 is obtained, similar to that obtained by Acker et al. (1992). The three pairs present 
enhanced [N n] emission, with mean [N ii]/Ha values of ~ 0.8 (A1-A2), ~ 0.15 (B1-B2), and ~ 
0.2 (C1-C2), which are a factor 2.5 to 13 higher than the mean ratio in the nebula. Furthermore, 
peak values of this ratio up to ~ 1.5 are found in A1-A2, and up to ~ 0.3 in B1-B2 and C1-C2. 
Common to all three pairs is that [N n]/Ha peaks at larger distances from the center than the 
[N n] emission peaks (see also Fig. 1). 

Table 1 lists different parameters of the bipolar pairs as deduced from the [N II] images and 
long-slit spectra. In order to deduce systemic velocities, we have first considered the centroid of 
each bipolar pair as representative of its own systemic velocity. The systemic velocity obtained 
for each pair is listed in Table 1. Secondly, we have obtained the centroid velocity of the [N n] 
emission line as representative of the systemic velocity of the whole nebula. In order to avoid the 
contribution of the bipolar pairs, the centroid velocity was obtained from the middle point of the 
spatially integrated emission at 10% level of the intensity peak. We obtain V S3/S (centroid) ~ 168 
kms -1 and note that an identical value is obtained if the middle point is considered at 5% or 15% 
levels of the intensity peak. The values for V SJ/S (centroid) , V sys (Bl-B2) and V sys (Cl-C2) agree 
within the absolute errors with the value quoted by Bianchi (1992) and Durand et al. (1998) for 
the systemic velocity of IC4846. The systemic velocity of B1-B2 and C1-C2 are almost identical 
to each other; both are slightly blueshifted (~ 3.5 kms -1 ) with respect to the centroid velocity. 
In the case of A1-A2, its systemic velocity is noticeably blueshifted by ~ 10 kms -1 with respect 
to the centroid velocity. We will discuss in § 4 possible interpretations for these systemic velocity 
shifts. In the following, radial velocities of each pair will be quoted with respect to its own systemic 
velocity. 

3.2. Properties of the collimated components 

IC4846 displays three pairs of structures, A1-A2, B1-B2, and C1-C2, which are highly evoca- 
tive of collimated outflows. We describe as follows the properties of these components. 

A1-A2 appear detached and off-axis (A(PA) ~ 43°) with respect to the major axis of the 
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inner ellipse. Their orientation changes by ~ 5°-10°. A1-A2 can be traced along ~ 6 arcsec, 
which represents a considerable fraction of the total nebular size. An upper limit of ~ 1 arcsec is 
obtained for their lateral size. The radial velocity systematically varies along A1-A2: a minimum 
velocity of ~ ± 36 kms -1 is observed at the minimum angular distance; the velocity then steadily 
increases with the distance up to ~ ± 52 kms^ 1 ; finally, it slightly decreases at the tips of the 
features. The velocity width is small as compared with the velocity difference between Al and A2 
(see Table 1), which implies a high degree of collimation. According to their properties, A1-A2 are 
highly suggestive of a precessing bipolar jet. 

B1-B2 present a rather moderate radial velocity and the velocity width is small and comparable 
to that of A1-A2 (Table 1). The small velocity width and the compactness indicate that B1-B2 
are highly collimated outflows. 

The radial velocity of C1-C2 does not depend on the PA (Fig. 3, Table 1), although some 
small variations (~ 7 kms -1 ) within C1-C2 are observed. Possible interpretations of this pair will 
be discussed in § 4. 



4. Discussion 

The observations of IC4846 have disclosed a complex structure composed of three different 
elements: (1) two apparent elongated shells with different orientations, (2) an attached circular 
shell and (3) three pairs of point-symmetric components at different orientations with typical prop- 
erties of collimated outflows. In the following we discuss about the nature, models and possible 
interpretations of these components. 

As already mentioned, there is no kinematic evidence for the elliptical shells in the PV maps 
of the [N n] line, although the extended emission observed in the PV maps may be related to 
these shells. We also note that the major axis of the outer ellipse is only defined by A1-A2. 
In fact, without A1-A2 the outer ellipse could not be recognized. Images at higher resolution 
would be necessary to elucidate the precise geometry of these shells. While A1-A2 and B1-B2 can 
be interpreted as precessing bipolar outflows, the interpretation of the pair CI-C2 is somewhat 
difficult. Its kinematic properties and the elongation perpendicular to the radial direction suggest 
that C1-C2 could represent collimated ejections in which the orientation of the ejection axis has 
remarkably changed over the ejection time. However, we cannot completely rule out that C1-C2 
represents a ring-like structure containing two bright opposite regions. The attached circular shell 
probably represents a previous slow wind ejected by the AGB progenitor of IC4846. It is similar 
to other attached shells observed in PNe and its formation may be related to the dynamics of 
ionization fronts (see Mellema 1995 and references therein). 

In the scenario that AI-A2 are precessing bipolar jets, we can obtain estimates for the inclina- 
tion angle and velocity by assuming that the expansion velocity is constant and that the variation 
of the radial velocity is due to a systematic change of the angle between the observer and the 
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precessing jet along a precession cone with aperture angle ~ 5°-10°. The observed minimum and 
maximum radial velocities are assumed to be related to the far and rear side of the precession cone, 
respectively. Using the model for a biconical flow described by Miranda & Solf (1992), we obtain 
an inclination angle for the precession axis of ~ 13°-26° with respect to the plane of the sky, and 
a velocity of ~ 100-190 kms^ 1 . The (mean) kinematic age is ~ 200xD(kpc) yr, where D is the 
distance to IC4846, which is extremely uncertain with published values between 2 and 11 kpc (see 
Acker et al. 1992). The images also show that A1-A2 protrude from the circular shell. This result, 
the increase of the [N n] /Ha ratio at larger distances and the decrease of the radial velocity at the 
tips of A1-A2 suggest an interaction between the jet fronts and the circular shell, in a similar way 
to that reported in, e.g., NGC6891 (Guerrero et al. 2000). 

It is worthwhile to emphasize the extraordinary similarity between IC 4846 and the well-known 
PN NGC 6543 (see Miranda & Solf 1992 and Reed et al. 1999). The following structural components 
are common to these two PNe (we will use for NGC 6543 the nomenclature of Miranda & Solf 1992): 

1. Basic elongated shell(s). In NGC 6543 two shells are observed, an outer bipolar shell with a 
bright ring between the lobes and an inner elliptical one. In IC 4846 two elongated shells are 
suggested by the images (but see above). The possible existence of an equatorial ring should 
be confirmed by new data. 

2. Bipolar knots at the tips of the shell. The knots FF' in NGC 6543 correspond to B1-B2 in 
IC 4846. 

3. Precessing bipolar jets. The bipolar jets JJ' in NGC 6543 correspond to A1-A2 in IC 4846. 
The similarity between JJ' and A1-A2 is extraordinary. Both appear off-axis with respect to 
the orientation of the basic shell and are curved, with a change of orientation of 5°-10°. The 
radial velocity presents an identical variation along the jets with a steady increase up to a 
maximum radial velocity and a small decrease at their tips. Finally, an expansion velocity of 
~ 100-200 kms" 1 is obtained in both cases. 

4. Point- symmetric components perpendicular to the radial direction. C1-C2 in IC 4846 could 
correspond to DD' in NGC 6543. 

There are some differences between both objects, but these are most likely attributable to 
different orientations, ages and distances. Whatever these differences, the conspicuous kinematic 
and structural similarities lead us to conclude that NGC 6543 and IC 4846 have shared a similar 
formation history. 

An interesting result in IC 4846 is the systemic velocity shifts among the different components 
and the centroid of the nebula. In the case of B1-B2 and C1-C2, the shifts are only slightly larger 
than the relative errors and we cannot conclude whether they are real. In the case of A1-A2, 
however, the shift is noticeable and we consider it real. Both the remarkable spatial and velocity 
symmetry observed in A1-A2 exclude that this shift may be ascribed to different ejection velocity, 
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changes (deceleration) in the velocity in one of the components of the pair, different ejection angle 
or a combination of the three. A more plausible explanation is that the systemic velocity of A1-A2 
is indeed intrinsically different from that of the other components and, in particular, from that of 
the emission centroid. Similar shifts in the systemic velocity of collimated outflows and shell have 
been recently reported in the PN Hu2-1 (Miranda et al. 2000). Furthermore, both PNe share 
the presence of two systems of collimated outflows at different orientations, one of them projected 
along the nebular major axis and one off-axis, suggesting a parallel formation history. Moreover, 
inspection of the kinematic data of NGC 6543 (see Miranda & Solf 1992) suggests that a difference 
of ~ 7 kms -1 could exist between the systemic velocity of JJ' and the systemic velocity of the 
nebula as defined by the rest of components. In the case of Hu2-1, the shifts of systemic velocity 
have been interpreted as a direct evidence of orbital motion of the ejection source in a binary central 
star (Miranda et al. 2000). The velocity shifts detected in IC4846 fit also into this framework. If 
A1-A2 are ejected from a star in a binary system, the space velocity of the A1-A2 will contain 
an additional velocity component due to the orbital motion, besides to the component due to the 
own ejection velocity. Accordingly, the systemic velocity of A1-A2 will be shifted from that of the 
emission centroid of the nebula (see Miranda et al. 2000 for details). We note that the same would 
apply to B1-B2 and C1-C2 if their observed velocity shifts were real. 

The systemic velocity shift of A1-A2 (10 kms -1 ) is very similar to that found in Hu2-1 (8 
kms -1 ) and in NGC 6543 (7 kms~'). Consequently, the orbital parameters of their binary central 
stars should be rather similar as well, with an orbital separation < 30 AU and an orbital period 
< 100 yr, thus linking the central stars of these nebulae to symbiotic stars. These results strongly 
suggest that the central star of IC4846, similarly to that of Hu2-1 and NGC 6543, is a binary 
star in which mass transfer from a red giant leads to the formation of an accretion disk about a 
secondary (Mastrodemos &; Morris 1998, 1999), from which collimated ejections may be generated. 

5. Conclusions 

Narrow-band imaging and high-resolution spectroscopy of the PN IC4846 in the [N n]A6583 
and Ha lines have revealed a complex structure similar in many respects to that of NGC 6543. 
The most remarkable components in IC 4846 are three pairs of point-symmetric structures oriented 
at different directions. Two (and probably the three) pairs represent collimated bipolar outflows. 
The kinematics and morphology of the pair observed at PA 54° (A1-A2) are highly suggestive of a 
precessing jet moving at ~ 100-200 kms -1 . Furthermore, the velocity shift between their systemic 
velocity and that of the nebula can be interpreted as a result of orbital motion of the ejection source 
in a binary central star. This kind of systemic velocity differences has also been observed in Hu2-1 
and could be present in NGC 6543. Thus, the formation of IC4846 (and Hu2-1 and NGC 6543) 
could be related to physical processes in the own evolution of interacting binary stars. 
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Tablel: Parameters of the [N n] bipolar pairs in IC 4846° 



Pair 


5\ r b 


AV C 


V sys (LSR) d 


5V sys e 


5Xf 






(kms -1 ) 


(kms -1 ) 


(km s _1 ) 


(kms -1 ) 


(arcsec) 


(degrees) 


A1-A2 


96 


20 


158 


-10 


7.4 


54 


B1-B2 


26 


18 


164 


-4 


5.4 


11 


C1-C2 


38 


22 


165 


-3 


1.8 


25 



a Obtained from the images and long-slit [N n] spectra shown in Figs. 1 and 3, respectively. 
b Separation in radial velocity between intensity peaks of the components of each pair. 
c Velocity dispersion (FWHM) corrected of instrumental resolution. 
d Systemic velocity. 

e Maximum difference between V sys (LSR) for each pair and V sys (LSR) for the centroid of the [N n] 
emission line =168 kms -1 

f Projected angular separation between the intensity peak of the components of each pair . 
9 Orientation (position angle). 
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Fig. 1. — Grey-scale maps of the [N n]A6583 (upper row), Ha (middle row) and [N n]A6583/Ha 
(lower row) images of IC4846. The [N n] and Ha images are presented in three different contrasts 
in order to show the faint and bright nebular regions. The observed slit positions at PAs 11°, 54° 
and 140° are indicated in one of the [N n] maps. In the [N n]/Ha map, black indicates stronger 
[N n] emission. North is up, east to the left. The size of the field shown is 14x14 arcsec 2 . 
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Fig. 2. — (upper row) Contour plots of the images of IC4846 presented in Fig. 1. The pairs of 
bright [N n] components are labeled A1-A2, B1-B2 and C1-C2. (lower row) Deconvolved Ha and 
[N n] images at a resolution of 0.5 arcsec. The contours are logarithmic separated by a factor 1.32 
in intensity. 
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Fig. 3. — (upper row) Contour plots of the long-slit [N n]A6583 spectra at three different PAs (upper 
left) in a position-velocity representation. The contours are logarithmic separated by a factor 1.41 
in intensity. The identified components are labeled (see Fig. 2). The horizontal line corresponds to 
the mean position of C1-C2. The vertical lines correspond to the systemic velocity of the different 
pairs (continuous lines) and to that of the emission centroid (dashed line) (see text), (lower row) 
Grey-scale maps of the [N n]A6583/Ha ratio deduced from the long-slit spectra of IC4846 at PAs 
11° and 54°. Black indicates stronger [N n] emission. 



